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ABSTRACT 



We study the plasma correlation effects on nonresonant thermonuclear reactions of carbon and oxygen in the interiors of white dwarfs 
and liquid envelopes of neutron stars. We examine the effects of electron screening on thermodynamic enhancement of thermonuclear 
reactions in dense plasmas beyond the linear mixing rule. Using these improved enhancement factors, we calculate carbon and oxygen 
ignition curves in white dwarfs and neutron stars. The energy balance and ignition conditions in neutron star envelopes are evaluated, 
taking their detailed thermal structure into account. The result is compared to the simplified "one-zone model," which is routinely 
used in the literature. We also consider the effect of strong magnetic fields on the ignition curves in the ocean of magnetars. 

Key words, dense matter - nuclear reactions, nucleosynthesis, abundances - stars: interiors - stars: neutron - white dwarfs 



1. Introduction 

Thermonuclear reactions play a crucial role in stellar evolu- 
tion scenarios. In particular, they are important in white dwarfs 
and neutron stars. Ignition of degenerate carbon and oxygen 
in the interiors of white dwa r fs gives rise to Type la super - 
novae (iHo vie & Fowled 1 1960t iHillebrandt & Niemevej I2000). 
Nuclear reactions in degenerate envelopes of accreting neutron 
stars are responsible for X-ray bursts and for the overall chemical 
and thermal structure of the envelopes dFushiki & Lambl[l987t 
Brown & Bildsten 1998). Nuclear fusion rates in the interiors 
of d egenerate stars can be significantly enhanced over the bi- 
nary iGamowl (119281) rates b ecause of the man y-body screening 
effect in the dense plasma (S chatzman 1948l_for reviews, see 
lYakovlev & ShalvbkovH9 89 and I chimarull993l) . 

The screening in degenerate matter is usually treated under 
the assumption that the electron gas can be considered as a uni- 
form "rigid" background. The influence of the electron polariza- 
tion on the enhancement of nuclear rea ction rates has been stud- 
ied in some detail in several papers dSalpeteiHl954l: lltoh et alJ 
1977t lYakovlev & Shalvbkovll 19891: ISahrling & Chabrietl ll99S 



Kitamura 2000), which confirms that the electron screening ef- 



fect is rather weak in degenerate matter. At the time of these 
studies, uncertainties in the reaction rates due to other factors, 
viz. quantum effects and deviations from the linear mixing rule 
(LMR) in strongly coupled plasmas, as well as theoretical uncer- 
tainties in the nuclear effective potentials at short distances, were 
more important than the polarizable electron-screening effects. 

The mentioned uncertainties have been substantially reduced 
in recent years. There has been significant progres s in treating 
cross sections of binary nuc lear fusion reactions dBeard et all 
1201 0l and references therein). Ya kovlev et al.l(l2010h constructed 



an analytic model for calculating these cross sections, which 
accurately describes the data and p arametrized it for a num- 
ber of C, O, Mg, and Ne isot opes. iPollock & Militzerl {2004) 
and iMilitzer & Pollock! d2005l) used the path-integral Monte 
Carlo (PIMC) method to determine contact probabilities of re- 
acting nuclei for one-component plasma (OCP) with emphasis 
on many-body quantum e ffects (these c a lculations supersede the 
previous PIMC study bv IOgatal[l997h . IChugunov etail d2007l) 
compared these PIMC results to semiclassical calculations and 
find good agreement between the two approaches at tempera- 
tures higher than about one fifth of the ion plasma tempera- 
ture. These authors also obtained a simple parametrization of 
the reaction rates wi t h allow ance for the ion quantum effects. 
IChugunov & DeWittl (12009a) extended these results to reactions 
between different nuclei and suggested an analytic expression 
for re action rates in multicompone nt ion mixtures, based on the 
LMR. Chuguno v"& DeWittl d2009b) used extensive Monte Carlo 
simulations and discuss corrections to the LMR for the plasma- 
screening function in strongly coupled binary ionic mixtures. 
They also propose an analytic formula for the screening func- 
tion in ion mixtures. 
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IChugunov & DeWittl (l2009albl) have employed the model of 
rigid electron background. In this paper we demonstrate that the 
electron screening effects are not negligible compared to the 
other improvements considered in recent publications. We de- 
rive a simple analytic formula for a quick evaluation of these 
effects. We also calculate the ignition curves for carbon, oxygen, 
and their mixtures. We consider plasma cooling by heat conduc- 
tion and different neutrino emission mechanisms, which evac- 
uate the heat released in nuclear burning, thereby determining 
the ignition curve. The account of the heat diffusion is taken by 
detailed calculation of the thermal structure of neutron star en- 
velopes and corresponding heat fluxes. The result is compared 
to the simplified "one-zone approximation" dBrown & Bildstenl 
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1998: lCumming & Bildster]l200ll:lGasques et al.ll2007l) . Finally, 
we consider the effects of strong magnetic fields on the ignition 
curves in neutron star envelopes. 

In Sect. [2] we compare different approximations for the en- 
hancement factors and study the effect of electron screening. 
In Sect. [3] we calculate carbon and oxygen ignition curves in 
degenerate stars, using the state-of-the-art treatment of carbon 
and oxygen fusion reactions, neutrino emission mechanisms, 
and heat conduction with allowance for strong magnetic fields. 
Results are summarized in Sect. [4] 



2. Enhancement factors 

2.1. Classical theory and modern approximations 

It is customary to write the cross section of binar y nuclear fusion 
reactions in the form (e.g., Yakovlev et al. 2010) 



oiE) 



■S(E), 



(1) 



where E is the center-of-mass kinetic energy of the reacting nu- 
clei "1" and "2", 



n = 



E R = 



{ZiZ 2 e l Ymi2 
2H 2 



(2) 



Zje is the charge of nucleus "/'", e is the elementary charge, 
= m\m.2l{m\ + 1112) is the reduced mass, and S(E) is a func- 
tion called "astrophysical factor." For the Boltzmann distribution 
of nuclei, the reaction rate (the number of fusion events per unit 
time in unit volume) in the absence of screening is given by 



R l2 = Wi2«i«2 



7TOT12 



T 3 



1/2 



r 

Jo 



e-^- E/T S(E)dE, 



(3) 



where rtj is the number density of the ions of type "/', T is tem- 
perature in energy units, and the factor w\2 accounts for statis- 
tics: w\2 = j, if nuclei "1" and "2" are identical; otherwise 
W12 = 1. If T is small, then the integrand in Eq. © is strongly 
peaked at the energy E^ _ (t^ErT 2 ) 1 ^ 3 , and the integral can be 
evaluated as (Salpeter & Van Hornll 1969b 



Jf e -W/T S (E)dE*\j-TE pk 



1/2 



S(£ pk )e" 



where 



r = 3(n 2 E R /T) 



1/3 



(4) 



(5) 



Approximation © is valid, if r » 1 . 

In order to take the plasma screening effects into account, 
it is convenient to write the radial pair-distribution function for 
ions in the form 



i 2 (r) = exp - 



ZiZ 2 e" 



rT 



exp^— 



(6) 



where the first factor is the Boltzmann formula for an ideal gas, 
while the second one shows how the probability of separation of 
two chosen ions is affected by the surrounding plasma particles. 
The functio n H\2{r) is often ca lled screening potential of the 
plasma (e.g jDeWitt et al.lll973h . 

Along with the customary ion sphere radii aj = 

11/3 



Coulomb coupling parameters Fj = (Zje) 2 /ajT, it is convenient 
to introduce parameters 



r,, = 



Z x Z 2 e 2 



an = 



a\ + a.2 



(7) 



Provided that H\2(r) varies slowly on the scale of the classi- 
cal turning po int distance, which requires that 3Fi2/r <k 1 
(Ilchimarulll993h . the screened reaction rate is appro ximate! 
given by R\z& h , where the enhancement exponent is (Salpefc 
fl954h 

h = H 12 (0)/r, (8) 



ter 



and R\2 is giv en by Eq. ® with replace ment of S(E) by 
S(E + gi 2 (0 )) (IChugunov & DeWitt) l2009i. As discussed by 
iMitlerl d 19771) and lltoh et alj I d 19771) . approximation <(8j needs to 
be corrected at higher densities, where 3>T 12 It is not small; in the 
latter case, the quantum effects on ion motion become significant 
dJancovicilll977l:lAlastuev & Jancovicllll978l) . 

The Helmholtz free energy F(V, T; {Nj};N e ) depends on the 
numbers Nj = njV of ions of all kinds, the number of electrons 
N e = n e V, volume V and temperature T. We write it in the form 
F = Fa + F ex , where F^ is the free energy of the ensemble of 
noninteracting ions and electrons, and F ex is the excess free en- 
ergy that accounts for the interactions. In this paper we consider 
only n eutral plasmas, so t hat n e = TjJijZ j. One can rigorously 
prove dDeWitt et al.ll973Lljan"covicifl977l) that H 12 (Q) equals the 
difference between the excess free energies before and after an 
individual act of fusion. In the thermodyna mic limit this gives 
the relation (cf . Ilchimaru & K itamura 1996) 



(d d d \ 



(9) 



where n mn = 2; «> is the total number density of ions, including 
number density «3 of composite nuclei, which have charge num- 
ber Z3 = Z\ + Z2 and mass OT3 » m\ + m.2, and / ex = F ex /rii on VT 
is the normalized excess energy. 

In strongly coupled Coulomb plasma mixtures of clas- 
sical ions and degenerate e l ectrons, the LMR i s fulf illed 
(lHansen & Vieillefossd I1976L IChabrier & Ashcrofj Il990h . so 
that / ex ~ f\ m , where 



/im({«j}, n e , T) = 2^ xjfj(n e , T). 



(10) 



Here, xj = n.jjn\ m denotes the number fractions and fj(n e , T) is 
the normalized excess free energy / ex for a plasma containing 
only the jth type of ions. Accurate analytic expressions for fj 
in the Coulomb liquid have been derived in our previous work 
dPotekhin & Chabrierl2000h . and we use these expressions here- 
after. It follows from Eqs. (0 and (TTOb that the enhancement 
exponent in the LMR approximation is 



him = fi(n e , T) + f 2 (n e , T) - f 3 (n e , T). 



(11) 



In the approximation of rigid electron background, this reduces 
to 



Vii = /ii(ri) + /ii(r 2 )-/ii(r 3 ), 



(12) 



where /u(r) is the normalized excess free energy of the OCR 
In the io n sphere appro ximation, fa(T) = -0.9 T, hence h\ m ji 
becomes (Salp eterll 19541) 



(3Z;/47m e ) 1/J , where n e is the electron number density, and hs = 0.9 (r3 —F\ -I"^). 



(13) 



2 



A. Y. Potekhin and G. Chabrier: Thermonuclear fusion in dense stars 



We no te, in passing, that Eq. (169) of I Yakovlev & Shalvbkovl 
(1989) recovers this equation with a factor 1.055 instead of 0.9. 
We will see, however, that the factor 0.9 provides a much better 
approximation to the accurate screening function. 

The LMR is not exact, and it becomes progressively inaccu- 
rate with decreasing Tj. When T, <g 1 for all /', the e xcess free 
energy F ex is described by the iDebve & Hiickell (1 19231) approx- 
imation, Foh = -VT/l2nD 3 , where D is the screening length. 
For nondegenerate electrons and ions, 



D 



-2 



d: 2 = 



Ane 2 



n e , D. 



_ Ane 2 y 1 
_ T 2—i n 



z 2 



(14) 
(15) 



In the approximation of rigid electron background, D = D mn . 
In the Debye-Hiickel approx imation, Eq. (0 yields for the en- 
hancement exponent (Salpeteri l 19541) 



ZiZ? 



1DH 



DT 



(16) 



For an arbitrary degree of degeneracy (but at not too strong 
Coulomb coupling; see Ch abrierill990l) . the screened interac- 
tion between ions is approximately described by a Yukawa po- 
tential, Z\Zie 2 er r l D I V, with D = k JF for electron screening or 

D - (k 2 F + DT^j ' for electron and ion screening. Here, 



k TF = Ane dn e /d/J. e 



(17) 



is the Thomas-Fermi wave number, and is the chemi- 
cal potential of Fermi ga s of electrons. Using Eq. (24) of 
IChabrier & Potekhinl (fr998). one can write kjp(n e , T) in analytic 
form. We note, however, that the Yukawa model corresponds to 
the Thomas-Fermi limit, e(k) ~ 1 + (kjp/k) 2 , for the static dielec- 
tric f unction e(k), which ma y only be justified at k <k kjp (see, 
e.g., Galam & Hansen 1976). Therefore, this model is inappro- 
priate at short distances (i.e., large wavenumbers k). In particu- 
lar, it is not applicable for the evaluation of the screening poten- 
tial at zero separation, Hi ?(0). 

For the general case, Salpete r & Van Hornl d 19691) proposed 
the following interpolation between the Debye-Hiickel and 
strong-coupling limits: 



h 



S«DH 



SVH 



lb 2 + h 2 
\"s ^ "dh 



(18) 



where /is and /idh are given by Eqs. ( fT3l — <TT~6l> . 

Another analytic approximation for th e enhancement fac- 
tor be yond the LMR was constructed by Chug unov & DeWitH 
(2009b), based on Monte Carlo simulation results for the rigid 
background model. 

These analytic approximations can be compared to the result 
given exactly by Eq. We write the normalized excess free en- 
ergy in the form / ex = fi m + / mix , where fi m is given by Eq. ([T0j>, 
and f m ix({xj], {Zj}; n e , T) is the correctio n to the LMR, whic h 
was recently obtained in analytic form (iPotekhin et alJ 12009). 
Then, from Eq. (0, we obtain the enhancement exponent 



ho = h 



ho 



d/mix(*l + X 2 + X 3 - £) 



(19) 



f=0 



where h\ m is given by Eq. (fTTT i. Figure [T] shows enhancement 
factor s for 12 C fusion in differen t approximations, normalized 
to the Sal peter & Van Hornl (IT969) enhancement factor approxi- 
mation ( fT8l . Here, we intentionally neglect ion quantum effects 
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Fig. 1. Plasma enhancement exponents for carbon fusion re- 
actions in different approximations, neglecting the ion quan- 
tum effects, normalized to the enhancement exponent given by 
Eq. (list , as functions of mass density for 2 isotherms (top panel: 
T = 10 8 K; bottom panel: T = 10 9 K). Dot-dashed lines 
correspond to the electron rigid background model. Namely, 
dot-short-dashed lines show ha given by Eq. ( TT~9b for a rigid 
backg round, while dot-long-dashed lines show the approxima- 
tion of lChugunov & DeWitj (l2009bl) . The other lines correspond 
to the polarizable electron background case: dotted line: LMR 
[Eq. (fTTVl: solid line: h [Eq. (fl9Vl; dashed line: /i appr [Eq. (f2TT>1 . 



and postpone their discussion to Sect. 12.31 We compare the an- 
alytic expressions for the OC P (thus rigid background: the fit of 
IChugunov & DeWitlll2009bl and the result of using Eq. ([T9]l for 
a rigid bac kground, i.e. wit h h\ m replaced by h\ m ^ and /| n ; x given 
by the fit of lPotekhin etail (|2009) for the rigid background case. 
Hereafter, this approximation will be denoted ha. Comparison 
of the two dot-dashed curves shows that these approximations 
agree with each other within typically 2%. 



2.2. Electron screening 

In Fig. Q] we compare the enhancement factors obtained usin g 
Eq. (Q2]l, where / mix is given by the fit of lPotekhin et alJ (|2009), 
for the case of a polarizable electron background, and for the 
OCP. This comparison illustrates the contribution of electron gas 
polarization to the screening exponent under the present condi- 
tions. Additionally we show him, given by Eq. (fTTT i. with the elec- 
tron polarization taken into account in fj(n e , T). We see that the 
correct enhancement factor differs appreciably from the LMR re- 
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suit in the low-density regime, and the difference increases with 
temperature. 

The electron screening contribution can be quickly estimated 
as follows. We define the effective ion screening length for the 
reacting nuclei as 



Ao„ = ^Dl n + (0.6 a 12 ) 2 . 



(20) 



At low densities, where the Debye-Hiickel theory is appropriate, 
Aon approaches D mn while at high densities it is proportional to 
ayi. The numerical factor 0.6 is the only fitting parameter. Then 
the approximation for ho reads as 



h- dpP r = huy[l+k 2 [F D* 0B . 



(21) 



The result is also illustrated in Fig.Q] 



2.3. Quantum effects 



|Jancovicil(ll977l) and lAlastuev & Jancovicil(ll978l) examined the 
short-range behavior of the internuclear correlation functions 
and showed that the quantum effects for the ions decrease the en- 
hancement factor, which we write as e h i . They developed a per- 
turbation expansion of the enhancement exponent /z q in powers 
of the para meter (3T/t). This th e ory is appli cable at T » 1 and 
3T/r < 1. iPollock & Militzeri (120041) and iMilitzer & Pollock! 
(2005) performed PIMC calculations of the contact probabil- 
ities in the quan tum regime. They confirm the conclusions of 
iJancovicil (1 19771) and Alastue v" & Jancovicl (Il978l) and extended 
numerical resu lts beyond the a p plicab ility range of the pertur- 
bation theory. IChugunov et alj ((2007) find that the results of 
Militzer & Poll ock! (12005^ agree with semiclassical calculations 
and suggest an analytic parametrization of the reaction rates 
that accounts for the quantum effects in an OCP. A similar 
parametriz ation for multicomponent mi xtures has recently been 
derived by IChugunov & DeW itt d2009al) . in the LMR approxi- 
mation. They find that the quantum effects can be described by 



the use of Eq. (fT2l i with the substitution^ Tj 
where 



1 1/3 



fi 2 = [l + c, (3r 12 /r) + c 2 (3r 12 /r) 2 + c 3 (3r 12 /r) 3 ]"~ , (22) 
ci = 0.013 y 2 , c 2 = 0.406 y 014 , c 3 = 0.062 y 019 + 1.8/r 



12, 



t and Fi2 are defined by Eqs. (0 and ©, and y — 4ZiZ 2 /(Zi + 
Z 2 ) 2 (seeFig.H. 

Expansions of the fitting functi ons of lChugunov et alj d2007l) 
and IChugunov & DeWittl (l2009a|) in Tavlor series do not re- 
cover the perturbation series of Alastuev & Jancovicll (Il978l) . 
This mismatch, however, is probably unimportant (unless one 
is interested in the second and higher derivatives of h), because 
the numerical agreement with the Alastuey-Jancovici results in 
their validity domain (3F/r < 1) is quite good. 

To include the quantum effects in the general case, we multi- 
ply the classical expression (e.g., Eqs. ( fT9l l or (Eq. (ET})) by the 
quantum decreasing factor q = /ji m , u 7/!i m> ii, where h\ m ^ is given 
by Eq. (TPS} , and h\ m \\ results from the replacement of Tj by Tj 
in Eq. (fT2l . The function hg ■ q is also plotted in Fig. [2] 



1 Here a typo in Chueunov & DeWitt ( 2009a) is corrected. We are 
grateful to A.I. Chugunov for drawing our attention to another typo in 
the first version of this our paper. 



1.1 




0.9 



0.8 



h (ions + electrons) ^ 

LM (ions+electrons) ^ 

1.M (OCP) \ 

CDW(a) \ 

h q (ions+electrons+quant.) \ 



I ' ' ' ' I ' i_i i_l 1 i_i ' I ' 1 1 1 l_i i_i lJ ii'il ii'il 



■=! 1.05 




Fig. 2. Same enhancement exponents as in Fig.Q] but with other 
approximations. Solid, dotted, and dot-dashed lines correspond 
to results while neglecting ionic quantum effects on the enhance- 
ment factor: the dotted line corresponds to the LMR for OCP 
[Eq. (fT2ll. the dot-dashed line demonstrates LMR for the polar- 
izable electron background [Eq. (II H I. and the solid line shows 
ho beyond the LMR approximation [Eq. (fT9bl. The dashed lines 
illustrate the impact o f ionic quantum effec t s: the l ong-dashed 
line presents the fit of lChugunov & DeWittl (l2009ah . while the 
short-dashed line is the approximation qho, which includes both 
the ionic and electronic screening contributions and takes both 
the quantum effects and the deviations from the LMR into ac- 
count (see text). 



2.4. Discussion 

Figures Q] and [2] demonstrate that electron screening al- 
ways increases the value of the enhancement factor. This re- 
sult is intuitively expected, because allowance for additional 
screening particles augments the overa ll effect. It qualita- 



tively a grees also w it h the fi ndings of IS alpeter & Van Hor 



( 1969) . iDeWitt et alJ (fl973h . lYakovlev ^S halvbkovNT98 



lord 



Sahrling & Chabrier Tll 19981). a nd Kitamura (2000). The opposite 
result was claimed bv lPollock & Militzeri (120041) . who found that 
electron screening "reduces the enhancement effect." This con- 
fusion arises from their use of a Yukawa potential to describe the 
electron screening. We pointed out in Sect. 12. 1 I that the Yukawa 
model becomes incorrect at short distances; in particular, it is 
incapable of determi n ing th e contact probabilities between fus- 
ing nuclei. Ilchimarul (1 19931) mentioned two opposite effects of 
electron screening: first, the binary repulsive potentials between 
reacting nuclei are reduced by electrons ("short-range effect"), 



4 



A. Y. Potekhin and G. Chabrier: Thermonuclear fusion in dense stars 



which increases //12(C)); second, the reduction of particle interac- 
tions by the screening affects the many -body correlation function 
in such a way that it decreases //12(C)) ("long-range effect"). In 
real electron-ion plasmas (without the Yukawa approximation) 
the first effect overpowers the second one. The Yukawa model 
grasps the second effect, but misses the first, dominant one. 

In the high-density domain (p > 10 8 g cm 4 ), where the fu- 
sion ignition in dense stars occurs, both corrections to h due to 
electron screening and due to the deviations from the LMR (see 
Fig- III) become quite small, of a few percent or less. Moreover, 
since these two corrections have opposite signs, they can bal- 
ance each other out. This is illustrated in Fig. [2] where the accu- 
rate result is the same as in Fig. [1] and we show also the result 
of application of the LMR for the cases of polarizable electron 
background according to Eq. (ITTb and the rigid background with 
LMR approximation according to Eq. Ill 2b . 

At lower density, the corrections to the enhancement expo- 
nent h due to electron polarization and deviation from the LMR 
become relatively large, but h itself becomes rather small, so that 
its effect on the reaction rate is not very significant. This domain 
concerns partially degenerate objects, such as low-mass stars or 
brown dwarfs, and may affect the nuclear production of light ele- 
ments, such as deuterium, lithium, beryllium, etc., which provide 
tracers for the mass and/ or age determination of these objects 
(IChabrier & Baraffefe OOO). These situations will be examined in 
a forthcoming paper. 

One can note that the simple approximation of 
ISalpeter & Van Hornl £1969), which we have chosen for 
normalization, performs surprisingly well: in the presently 
explored p—T domain of astrophysical interest, it provides 
an accuracy of the enhancement factor better than 10%, i.e., 
better than 0.04 for log m <t(E). Some later approximations 
dYakovlev & Shalvbko4l989YlItoh et al.lll990h are not as good. 
(They would fall outside the frames of Fig. [2j). 

Since the uncertainties in the n uclear part of the r eac- 
tion rates are still larger (see, e.g., Ag uilera et alj 2006), m 
practice it appears suffici ent to use the approximation of 
Salpe ter & Van Hornl ( 1 19691) with the quantum correction: 



h * /i SVH q, 



(23) 



where hsvn and q are given in Sects. I2TTI and 12.31 respectively. 
Including further improvements to h is currently just a ques- 
tion of completeness. However, as we have seen in Sect. 12.21 
the electron-polarization corrections are generally comparable 
to or even larger than other corrections discussed in literature. 
Therefore, the electron screening should be taken into account in 
every treatment of the thermodynamic enhancement factor that 
goes beyond approximation ( 1231 . 

3. Ignition curves 

3.1. Nuclear heating and neutrino cooling 

The ignition curve is the line in the p—T plane that determines 
the highest densities and temperatures at which exothermic nu- 
clear reactions in the plasma can be stable against thermal run- 
away. It is determined by the balance between nuclear energy 
generation rate and local heat losses. In this subsection we focus 
on the case where the heat losses are mainly caused by neutrino 
emission, which is appr opriate in white dwarfs, e.g., f or model- 
ing supernova la events dHillebrandt & Niem ever 200J). 

The thermonuclear fusion rate is given by Eq. (O, w here fo r 
S (E) we substitute the parametrization of Yakovl ev et al.l(l2010l) . 
The energy release power per unit volume equals RnQn, where 
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Fig. 3. Ignition curves for carbon (lower/left group of lines) 
and oxygen (upper/right group of lines) in different approxima- 
tions. Solid lines are obtained using Eq. ( fT9l and the quantum 
correction according to Sect. 12.31 long-dashed lines: Eq. ( fT2l 
and the quantum correction; dot-dashed lines: Eq. ([TBI without 
any correction. The lines are plotted heavy if the characteris- 
tic nuclear fuel burning time fb U m is shorter than 1 Myr; the 
lines are plotted thin if 1 Myr < fb U m < 14 Gyr. In the do- 
main where fb U m > 14 G yr the lines are dotte d. For the carbon 
case we also plot the fit ( Pote khin et al. 2003) to the results of 
ISahrling & Chabrier! d 19981) (the short-dashed line). 



<2i2 is the energy releas e in a single fusion event. We use the 
values of Q n given bv iFowler etail (1 19751) : 13.931 MeV for 
I2 C + I2 C, 16.754 MeV for 12 C + 16 0, and 16.541 MeV for 

16q + 16q react i ons 

At very high densities, pyc nonuclear fusion due to zero-point 
ion vibrations dCameronll959 ) becomes more important than the 
thermonuclear o ne. We calculate it usin g Eq. (33) and Table II 
(the first line) of I Yakovlev et al.l (120061) and add it to the ther- 
monuclear rate. The ignition curves portrayed in the figures, 
however, never enter the domain of pycnonuclear burning. 

The dominant mechanism of energy loss due to neutrino 
emission differs, depending on the physical conditions. For the 
present conditions of interest, the main mechanisms are neutrino 
bremsstrahlung by electron scattering off nuclei, plasmon de- 
cay, and electron-positron annihilation. The ener gy loss rates are 
given, respectively, by Eqs. (76), (38), and (22) of lYakovlev et al.l 
(120011) . 

Figure [3] illustrates the ignition curves for 12 C+ 12 C and 
i6q + i6q reac tj ons m different approximations. Here, the igni- 
tion curves are obtained with the enhancement factors given 
by Eq. ( fT9l with the quantum correction included according 
to Sect. 12.31 and are copmp ared with the approximat ions of 
IChugunov & DeWittl (l2009al) and lSalpeter & Van Hornl dl969h . 
In the low-temperature region, where the characteristic burning 
time fbum that is required to consume 63% of the nuclear fuel ex- 
ceeds the Universe age, the ignition curves lose any astrophys- 
ical sense, because the burning becomes unrealistically slow, 
and also because the poorly known quantum effects become too 
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strong (see the discussion in Gasques e t alj2005llYakovlev et al.l 

We see that the electron screening slightly shifts the igni- 
tion curves to lower densities, a consequence of the increased 
enhancement factor, as mentio ned previously. This r e sult ag rees 
with the previou s findings by ISahrling & Chabrier] (1 19981) and 
Kitamura (2000)Q In the figure we have a l so plo tted the car- 
bon ignition curve o f ISahrling & Chabrierl dl998l) as fitted by 
Poteldrinetalj{2003), which is close to our current result. The 
difference is mainl y caused by modern improvements in the as- 
troph ysical factor (lYakovlev et al.ll2010l) and neutrino reaction 
rates dYakovlev et al.l [2001b . included in the present treatment. 
In AppendixlAlwe present a fit to the current carbon and oxygen 
ignition curves in a wide density range. 



3.2. Conductive cooling 

Nuclear burning in neutron stars occurs in relatively thin en- 
velopes, whose geometri c depth is not more than a few percent of 
the stellar radius (e.g., iBrown & Bildsten] [l998. and references 
therein). Thermal conductivity in these envelopes is high and 
thermal relaxation time is short, so that cooling by heat diffusion 
can stabilize nuclear burning beyond the ignition limits consid- 

ere d in Sect. 13. II 

lHansen & Van Hornl d 1975b estimated stability of H and He 
burning in envelopes of neutron stars by comparison of the char- 
acteristic time scale of thermonuclear heating with characteris- 
tic time scales for removal of ene rgy from the shel l by ra diative 
and conductive thermal diffusion. iFushiki & Lambl d 19871) intro- 
duced a differential criterion for such estimates. They defined 
the boundary of thermal stability of the nuclear burning from the 
condition 



Se mc (y,T) Se coo \(y,T) 



ST 



ST 



(24) 



where y is the column depth of the burning material (carbon or 
oxygen in our case), ST is a variation in temperature, Se nuc is the 
respective variation in the nuclear energy release rate per unit 
mass (e nuc = 2?x2 Qu/p), and Se coo \ is the variation of the cooling 
rate e coo i = -dF r (y)/dy, where F, is the outward radial heat flux 
through unit surface. The cooling rate e coo i depends on the first 
and second derivatives of the temperature, therefore the condi- 
tion (l24l) is not local, but depends on the overall temperature and 
density structure of the envelope. 

Thermal relaxation of the neutron star crust occurs on the 
scale of a few tens of years, while the relaxation time of the outer 
envelopes is still shor ter dLattimer et al.l Il994t iGnedin et alJ 
et al.l20 10i). Therefore, the outer envelopes are con- 
sidered quasistationary for most of the astrophysical problems, 
and t heir thermal structure is c alculated at st ationary equilib- 
rium ( Gud mundsson et aTl l 1983b . In particular, Fushi ki & Lambl 
dl987 ) applied this approximation to the sta bility analysis of H 
and H e she lls of accreting neutron stars, and lBrown & Bildstenl 
(1998|) and ICumming & Bildstenl d2001l) applied it to the sta- 
bility analysis of carbon shells. The problem is that Eq. (l24l > 
depends not only on the equilibrium quasistationary tempera- 
ture profile T(y), but also on the profile of its variation ST(y). 
In equilibrium dF r (y)/dy — e v - e nuc + Tds/dt, where e v is the 
neutrino emission power per unit mass, s is the entropy per unit 
mass, and t is time. (Here for simplicity we neglect the General 



2 The carbon ignition curve of Kitamura (2000) substantially differs 
from our results, because it was calculated under the assumption of a 
fixed neutrino emission power. 
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Fig. 4. Carbon ignition curves in the ocean of a neutron star 
with mass M - 1.4 M and radius 12 km, according to different 
models: solid line corresponds to Eq. (1261 1. long-dashed line to 
Eq. d25l ). short-dashed and dotted lines, respectively, to the orig - 
inal and modified approximations of Fushiki & Lambl (Q987), 
Eq. (1241 . For comparison, the ignition curve without account of 
conductive cooling is plotted by the dot-dashed line. 



Relativity corrections; for the c omplete set of accurate equations 
see, e.g., Richardson et ai ll979h In the quasistationary approxi- 
mation, one neglects ds/dt. As a result, were the variations ST(y) 
in equilibrium, one would have e coo i = e nuc - e v , which ex- 
hausts Eq. (12411 . Meanwhile, different nonequilibrium forms of 
ST(y) lead to different values of fen^i . IFushiki & Lambl dl987) 
assumed the global temperature variation, dST/dy = 0. For ex- 
ample, an assumption of the overall variation in the internal en- 
ergy would result in d<S77dy = -(ST/cp) dcp/dy, where cp(y) 
is the heat capacity per unit mass at constant pressure, taken at 
equilibrium. 

In modern l iteratu r e dBrown & Bildstenl Q998; 
Cumming & Bildstenl 1200 U iGasques et al.l 12007b . a simpli- 
fied "one-zone approximation" is used for the cooling rate in 
Eq. (|24|): 



ecooi =PKT/y 2 , 



(25) 



where k is thermal conductivity, y - P/g is the column depth 
in the plane-parallel nonrelativistic approximation, P is pres- 
sure, and g is the surface gravity. Definition d25l l is local and, 
therefore, free of uncertainties. Because it is local, the varia- 
tional derivatives Se/ST in Eq. ( l24l > can be replaced by the par- 
tial derivatives de/dT. However, it misses the information about 
the real dependence of conductivity on column depth, which is 
fraught with risk being inaccurate. 

We define the ignition curve in the outer envelope of neu- 
tron stars in the quasistationary approximation, taking into ac- 
count the detailed thermal and mechanical structure of the en- 
velope and assuming equilibrium variations of the temperature 
profile ST(y). For definiteness, let us suppose that ST > 0. For 
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a surface element dS of a shell confined between yi - y and 
yi - y + dy, an increase in temperature 5T leads to an increase 
in the nuclear energy release power 5e nuc dydS , an increase in 
neutrino energy emission power 6e v dydS , and changes of the 
heat flux through the outer and inner boundaries of the shell, 
SFi ^dS = 6F r (yi 2)dS . If 6e nnc ^ 6e v , the increase in the nu- 
clear heat release is compensated for by the increase in the en- 
ergy emission, and the nuclear burning is stable. 

We now consider the case where Se nuc > 6e v . If the sta- 
tionary equilibrium requires an increase in the heat income rate 
F2 dS through the inner boundary, i.e., 6F2 > 0, it means that 
the increase in the net heat release in the considered volume, 
6e tot = <Se nuc - 5e v , is overcompensated for by the increase in 
the outward flux to the stellar surface. In this case, the nuclear 
burning is stable. In the opposite case, where 8F2 < 0, 6e tot is 
not compensated by the conductive energy escape to the surface, 
so that thermal runaway occurs. Therefore, the largest column 
depth y, at which the burning can be stable, corresponds to the 
condition 

8F r (y)/8T(y) = 0. (26) 

At the column depth where the condition d26*b is satisfied, the net 
energy release 8e to i is balanced exactly by the increase in surface 
luminosity. Beyond this stationary point, thermal runaway starts. 

Figure [4] illustrates different approximations for the carbon 
ignition curve in the ocean of a typical neutron star with mass 
M - 1.4 M and radius 12 km. The carbon ignition curve with- 
out conductive cooling is calculated according to Eq. ( 1A.U . The 
other lines have been obtained by calculating series of tempera- 
ture profiles for carbon envelopes of the star, assuming different 
surface luminosities and applying different ignition conditions. 
The t hermal structure has b een calculated using the same code 
as in [ Kaminke retal] d2009), but the artificial heating model of 
Kaminkeret afl has been replaced by the nuclea r heating. The 
line corresponding to the iFushiki & Lambl dl987l) model shows 
a jump at a certain temperature, which corresponds to a switch 
from conductive to neutrino cooling. The switch signifies that 
at higher temperatures the constant variation 6T cannot pro- 
vide a powerful enough off-equilibrium cooling rate <Se coo i to 
compete with the neutrino cooling. Another functional choice 
of 8T(y), corresponding to a constant variation of internal en- 
ergy per mass, is shown by the dotted line. It displays a similar 
switch at a higher temperature, therefore the conductive cooling 
provided by this variation appears more efficient. The underes- 
timation of the ignition densities and temperatures by the two 
artificial forms of 5T(y) above the jump temperatures indicates 
that these functional variations are unstable on the characteristic 
cooling timescale. 

The line corresponding to the one-zone approximation d25] > 
is compared to the positions of the stationary points d26b . The 
latter two models qualitatively agree with each other, with the 
quantitative difference caused by the local approximation em- 
bedded in the one-zone model. 

A similar comparison of the models for the oxygen ignition 
curve gives similar results. For typical accreting neutron stars 
with low magnetic fields and effective temperatures from one 
to several MK, the ignition curves obtained in our calculations 
using the "stationary point" condition d26b can be reproduced by 
a modified one-zone approximation, where the right-hand side 
of Eq. d25l l is multiplied by a constant factor a (or equivalently, 
an effective thermal conductivity = ok is substituted). As the 
surface gravity varies from g » 0.7 x 10 14 cm s~ 2 (e.g., for stellar 
mass M = M & and radius about 15 km) to g m 4 x 10 14 cm s~ 2 
(e.g., for M — 2 M and radius of 10 km), the correction factor 
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Fig. 5. Carbon (lower/left group of lines) and oxygen (up- 
per/right group of lines) ignition curves in the ocean of a mag- 
netar with mass M - 1.4 M , radius 12 km, and magnetic field 
B = 10 14 G for different angles between the field lines and the 
normal to the surface: 6 B = 0°, 60°, 75°, and 85° (respectively, 
4 solid lines from top to bottom in each of the two groups). For 
comparison, the ignition curves given by Eq. ( IA.U are plotted by 
the dot-dashed lines, and the ignition curves without accounting 
for the magnetic field effects are shown by dashed lines. 

a varies from 0.18 to 0.16 in the case of carbon shell and from 
0.33 to 0.28 in the case of oxygen shell burning. 

3.3. The effect of strong magnetic field 

A strong magnetic field can affect conductivities and make the 
heat transport anisotropic, if the Hall parameter (the product of 
electron gyrofr equency and effective r elaxation time) is larger 
than one (e.g., Urp in & Yakovlevi ri980, and references therein). 
Therefore, it can affect the conductive cooling rate and shift the 
ignition curves in neutron star envelopes. We have evaluated the 
magnitude of this effect by calculating the temperature profiles 
for different magnetic field strengths B and inclinations 6b at the 
neutron star surface. 

In strongly magnetized neutron star envelopes, synchrotron 
neutrino emission becomes an i mportant energy sink. We calcu- 
late it according to Eq. (56) of lYakovlev et al.l d200ll) . Figure [5] 
illustrates a few examples. Here, a neutron star with M — 1 .4 M 
and radius 12 km is suppose d to possess a ma gnetic field B = 
10 14 G, typical of magnetars (Mereghet til2008l) . The lower lines 
are drawn for the carbon envelope, and the upper lines for the 
case of oxygen. In each of these groups of lines, the four solid 
curves (from top to bottom) are plotted for angles 6b between the 
field lines and the normal to the surface equal respectively to 0°, 
60°, 75°, and 85°. We do not show the case of a strictly tangen- 
tial field, because, for any reasonable distribution of the mag- 
netic field over the stellar surface, the heat transport becomes 
essentially two-dimensional around the surface spots where 6b 
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is close to 90°, which invalidates the employed approximation 
of radial (one-dimensional) transport. For comparison, we show 
the nonmagnetic ignition curves and the fusion-vs-neutrino bal- 
ance curves discussed in Sect. 13.11 (presumably appropriate for 
white dwarfs). 

When the field lines are nearly tangential to the surface, 
the strong magnetic field increases the conductive opacity, be- 
cause the heat transport perpendicular to the field is strongly sup- 
pressed. Such a field lowers the ignition curves and thus reduces 
the stability region. 

In the opposite case of field lines perpendicular to the sur- 
face, the strong magnetic field somewhat raises the ignition 
curves and expands the stability region. The latter result could 
not be obtained in the one-zone approximation ( l25l l. which re- 
lies on the local values of thermal conductivity, k, and ther- 
modynamic quantities p, T, and P. Indeed, in the p—T do- 
main shown in Fig. [5] thermal conductivity along the magnetic 
field is the same as without the field, because the field is non- 
quantizing at p » p B , where, for carbon-o xygen composi tions, 
p B = 1.41 x 10 7 CB/10 14 G) 3/2 g cm 3 (e.g.. lPotekhinll999t) . The 
equation of state also does not depend on B, if p » ps- Finally, 
the synchrotron neutrino emission cannot produce the observed 
effect, because the neutrino synchrotron power is less than 5% 
of thermonuclear power along the full ignition lines for the con- 
ditions appropriate to Fig. [5] The effect is caused by the opacity 
decrease in the heat-blanketing envelope at densities p < pg, 
where the field is strongly quantizing. The opacity reduction in 
this outer region makes the entire envelope more transparent and 
facilitates the heat escape to the surface. Thus, the effect is in- 
trinsically nonlocal. 

4. Conclusions 

We have studied the effects of electron screening on thermonu- 
clear reactions in dense plasmas and compared different approx- 
imations to determine plasma enhancement factors for the nu- 
clear fusion rates. The electron screening always increases the 
enhancement effect. The opposite conclusion, sometimes en- 
countered in the literature, comes from using the Yukawa po- 
tential model, which is inappropriate to calculating the contact 
probability for fusing nuclei. The electron screening correction, 
which we calculate using accurate analytic expressions for the 
free energy, can be satisfactorily described by a simple approxi- 
mation (ETT l. This correction, although small for the case of dense 
stars, has the same order of magnitude as other recently sug- 
gested corrections to the enhancement factor, so it needs to be 
properly calculated. 

Using our analytic formulae for the enhancement factors 
and state-of-the-art astrophysical factors for thermonuclear cross 
sections, we calculated carbon and oxygen ignition curves in de- 
generate stars. The ignition conditions in neutron star envelopes 
were evaluated, taking their detailed thermal structure into ac- 
count. Comparison of the results to customary simplified models 
demonstrates the restricted applicability of these latter. 

We also studied the effects of strong magnetic fields on the 
ignition curves in neutron star envelopes. These results show 
that the ignition surface shifts to lower densities in the stellar re- 
gions where the magnetic field is strongly inclined and to slightly 
higher densities in the regions of nearly radial magnetic field. 
The latter effect could not be obtained in the simplified one-zone 
model. For a magnetar, the shift in the ignition curve can be sim- 
ilar to the difference between the accurate and one-zone calcula- 
tions and larger than the correction due to the electron screening 
in the plasma. 



Table A.l. Parameters of Eq. ( 1A.U for the cases I6 0+ 16 reac- 
tion in pure oxygen plasma and 12 C+ 12 C reaction in a mixture 
of carbon and oxygen with number fractions xc and xo, respec- 
tively (0 x < 0.99, x c = 1 - xo). 



Parameter 



Pure 16 12 C+ I2 C in 12 C/ 16 mixture 



log 


ioPi.m (g cm" 3 ) 


10.6902 


9.76 + 0.025j o -0.47 In x c 


log 


io Pi (g cm ~ 3 ) 


7.2 


4.82 + 0.6x o -0.2 In x c 


log 


io P2 (g cm" 3 ) 


6.6 


2.8 + 2.2x 


log 


10*1 (K) 


9.161 


8.75 + 0.015xo - 0.033 In x c 


<>'() 




3.55 


4.2 






1.4 


0.5 + 3x 


"'2 




0.155 


0.084/ff! -0.0029 In x c 


a 3 




0.085 


-0.0053 In x c 



Acknowledgements. A.Y.P acknowledges useful discussions with Dima 
Yakovlev and partial support from the RFBR Grant 1 1-02-00253-a and the 
Russian Leading Scientific Schools program (Grant NSh-3769.2010.2). 



Appendix A: Fit for C and O ignition curves 

White dwarfs contain a mixture of carbon and oxygen nuclei. 
Therefore, it is of practical interest to determine ignition curves 
for such mixtures with different number fractions of carbon (xc) 
and oxygen (xo = 1 - xc). We have derived a fitting formula 
for ignition temperature as a function of mass density, which is 
relevant for densities p > 100 g ctrT 3 , temperatures T > 10 8 K, 
and oxygen number fractions xq ^ 0.99: 



7fit(p < Plim) 



Ti[i + (pi/prr 



{1 + [ao/ln(pi im /p)] 2 - 7 } - 2 [1 + (p 2 /p) 2 ] 



\2-\at, 



(A.l) 



where pii m , pi, pz, T\, and czq — as are fitting parameters, given 
in Table IaTTI 

At high p, the ignition curves determine the critical density 
of the ignition, rather than the critical temperature, as at smaller 
p. Therefore, there is no sense to measure the fit error by dif- 
ferences between the model (T&t) an d data (Tdat) values of tem- 
perature at a fixed p. Instead, we measure the fit error by the 
geometric distance in the logp-log(r 4 ) plane between numer- 
ical points (pdati Tfa t ) and the line Tt(p) (the fourth power of 
T is relevant because it is proportional to the luminosity). The 
fractional error is defined as 



e(Pdat) 



min 

p 



( T 4 (o) - T 4 ^ 2 

1 & 1 dat 

1 dat 



P - Pdat 
Pdat 



1/2 



(A.2) 



The maximum error max Pdat e(pdat) varies from 0.06 for xq - 0.5 
to 0.09 for xo = and xq - 0.99. For oxygen, maXp dil e(pd a t) = 
0.05 under the condition that T > 1.7 x 10 s K (at lower temper- 
atures the ignition curve is unreliable anyway). This is a good 
accuracy, given that the considered p and T 4 values span several 
orders of magnitude. 
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